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We study the statistical mechanics of multi-index matching problems where the quenched disorder 
is a geometric site disorder rather than a link disorder. A recently developed functional formalism 
is exploited which yields exact results in the finite temperature thermodynamic limit. Particular 
, attention is paid to the zero temperature limit of maximal matching problems where the method 

allows us to obtain the average value of the optimal match and also sheds light on the algorithmic 
heuristics leading to that optimal match. 



I. INTRODUCTION 



Matching problems have been studied both from the perspectives of combinatorial optimization 1] and statistical 
physics 01- In the statistical physics approach, interest lies in computing average properties of the matching path in 
stochastic versions of the problem. The most general form of the m-partite multi-index problem is stated as follows. 
Q • One has m sets of N objects and a matching consists of combining one element from each set into a molecule or match 
^ ' with m components. If in the i-th match the element we choose from set j, 1 < j < ?n, is denoted by ij, then the cost 
of the i-th match is given by the entry Ci-^...i^ of the cost matrix. The total cost of the matching is then given by the 
sum of the costs of the N individual matches. The underlying maximum/minimum optimization problem consists of 
^ finding the global match which maximizes/minimizes this cost. The bipartite matching problem with m — 2 has been 
most studied and some analytic progress has been made for so called "random link" models where the elements of 
] the cost matrix C are chosen independently from some probability distribution H, Q . 

Often matching problems have a geometrical origin as the points to be matched lie in some d-dimensional domain 
T). Euclidean matching problems are of practical importance, finding applications in air traffic control 5] and the 
manufacture of printed circuits (y\. The three dimensional Euclidean matching problem, both minimal and maximal 
has been studied numerically via Monte Carlo simulations and in particular simulated annealing J3| . In these geometric 
, problems the matching matrix C is determined from the location of the points using some metric on the space. As 
■ a result matrix elements are no longer independent (for example triangle inequalities exist for the case when the 
matching matrix depends on the Euclidean distance between the points). The random link versions can be regarded 
as approximations to Euclidean models, and in minimal bipartite matching this approximation seems rather good @ 
' O ' and corrections to account for triangular correlations have been computed Q. 

C , Recently, we have developed new techniques for analyzing Euclidean traveling salesman type problems (TSP) 
based on functional integration and a functional order parameter. This formalism does not yield the average value 
of the optimal path for the minimal TSP because it does not apply when the energy is rescaled to make the ground 
state energy extensive. However in the case of maximal TSP type problems the ground state energy is extensive 
without any energy rescaling and the average value of the optimal solution can be computed exactly. The method 
5^ ■ does not rely on the replica trick, however the results obtained show that for the Euclidean TSP the system is 
d ' replica symmetric, as is the case for random link calculations 0|- In this paper we show how the the formalism of 
[3 can be applied to multi-index matching problems at finite temperature and also how it can be used to obtain the 
average value of the optimal match in the maximal case. In addition the method gives us information concerning the 
algorithmic heuristic leading to the optimal path. The mathematical structure of the multi-index matching problem 
in our thermodynamic formalism is more complicated than that of the Euclidean TSP in that a functional order 
parameter appears for each of the m sets of points. However in the cases studied here the functional order parameters 
are the same for each set and we find no evidence of symmetry breaking. As a consequence of this symmetry we find 
that the thermodynamics of the bipartite matching problem and the Euclidean TSP are essentially equivalent. For 
certain other multi-index problems, with chain-like matching functions, we find relations between the average path 
lengths computed in these models and the TSP. 

The statistical mechanical cavity approach has recently been applied to minimal multi-index matching problems 
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in the context of independent links In this problem the energy (or equivalently the temperature) is scaled to 
ensure that the ground state energy, corresponding to the maximal match, is extensive. In contrast to the bipartite 
matching problem which has a replica symmetric solution, it was shown that the multi-index matching problem for 
more than two indices has a low temperature glassy phase characterized by replica symmetry breaking. It is therefore 
interesting to sec if the formalism developed in 9J which is exact for maximal optimization problems shows similar 
behavior when applied to multi-index matching. However in all of the cases studied here the system exhibits no such 
phase transition, suggesting that they appear to be non-glassy from the statistical mechanics point of view. 

II. MODEL AND GENERAL SOLUTION 

Here we define the general multi-index matching problem with site disorder. 

Consider the following m-partite matching problem with m — K + 1 sets S^°\ S^^\ • • • , S'^-^^ each set S*-"^ consists 
of N points {r^"-* , ■ • ■ , r j^' } distributed in a common domain I? in a space of dimension d. We assume that the points 
are independently distributed within the domain and that those in set S"*^"^ are distributed according to a probability 
density pa- An individual match consists of K + I points where one point in each set is matched with one point in 
each other set. Each match has an energy function associated with it which we denote by V{r'-°\r'^^\ ■ ■ ■ r*^^^). The 
multi- index optimization problem involves making N individual matches (Fig. Q), where each point is associated 
with one and only one matching, and optimizing the total match energy. The case oi K = 1 corresponds to the well 
known bipartite matching problem. A micro-state of the model is specified by {(i,a^^^ ■ ■ •Cj-^'')} where the point i 
of the set S^'^'^ having position rl"-* (1 < « < N)) is matched with r*-^;!, • • • ,r^^) in the other sets and CTj-"'' denotes 

the label of the element in set a chosen to be matched with the element i in set 0. The a^""^ are permutations on 
the N elements of 5'*^°'. The phase space for multi-index matching is thus the Cartesian product (Eat)^, where Ejv 
is the permutation group on N elements. The size of this phase space is (A^!)^ and the entropy is consequently 
super-extensive. 

The Hamiltonian of the system is given by adding up the energy for each individual match 

N 

H[a^,a2, • • • , ax] = ^ V{rf,rl% , • • • , r^fi,). (1) 
1=1 

and the partition function for sets of N points is given by 

Zn= Y1 e^p{-l3H[ai,a2,---aK]), (2) 

cr(l),cr(2),...cr(iC) 

where /3 = 1 /T and T is the canonical temperature of the system. 

The form of V can be arbitrary, but a symmetric potential function where V{ro, • • • , tk) is left invariant by any 
rearrangement of its arguments is natural in this context. A simple way of constructing such a potential is to consider 
the potential to be the, totally connected, sum of pairwise matching potentials between all pairs as shown in Fig. |(5J: 

K 

V{ro,ri---rK)^ ^ Vpira,n). (3) 

a=l,6=l 

However, to make analytic progress, much of our analysis considers potentials that are pairwise but are not sym- 
metric. These potentials allow the match to be represented as an ordered path through points from each successive 
set as illustrated in Fig. (0). 

K 

F(ro,ri...rK) =^Fp(r,_i,rJ. (4) 

a=l 

Besides this open model, a generalization is to add Vp(rx,ro), to arrive at the situation shown in Fig. I^J and we 
shall term this a closed or cyclic model. 

Because we emphasize the geometric interpretation of pairwise matching, we will frequently take the pair potential 
Vp to be the Euclidean distance Vp(r,_i, rj) = jr^-i — r^j, though we also treat the quadratic form which seems more 
amenable to analytic approaches. The case K = 1 corresponds to the Euclidean bipartite matching problem. The 
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FIG. 1: Full match for three sets (white, black, gray) of three points (m = 3, A*' = 3). 




FIG. 2: Individual match between four points (each in a different set white, gray, black and stripes), showing a fully symmetric 
and fully connected cost function. 

cyclic, K = 2 tripartite case is symmetric and corresponds to the triangle interaction which amounts to the total 
length needed to form a triangle between three points. Another symmetric possibility for K = 2, tripartite matching, 
that is not of the pairwise form, is the star interaction where the three points are connected by three lines emanating 
from the point which would be their center of mass (for equal masses) 

Vstar{Ya,Yi,Y2) = |ro - R| + |ri - R| + \y2 - R|, (5) 

where 

R= ^(ro + ri+r2). (6) 

As pointed out in the disorder in this type of problem can be encoded in the unaveraged density of the quenched 
points in the domain T) 

1 ^ 

1=1 

note that the value of pa averaged over the disorder is by definition pa- 




FIG. 3: Individual match between four points (each in a different set white, gray, black and stripes), showing a open chain cost 
function. 
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FIG. 4: Individual match between four points (each in a different set white, gray, black and stripes), showing a closed chain 
cost function. 



Imposing this constraint we may discard the sum over permutations and write 



iVp„(r)-^<5(rf)-r) 



i=l 



exp -/35:y(rf),...,rf)) . 



(8) 



In the above formulation without the delta function constraints the matches chosen can use any point in the domain 
V, however the delta function constraint allows the use of only the points which are available in the sets Sa- The 
delta function constraints which are present at each point r in the domain T> are now expressed as a functional Fourier 
integral (the /U integration is along the imaginary axis) 

Zn = C'^J l[d[f,,]l[dr\''h^p(N J cir^M,(r)pJr)]exp(-^M,(r(''))-/3E^(^i"V- - 'rl""^)) • (9) 

a i,a \ a / \ ai i / 

The integrals over the dynamical variables r^"^ which are taken to be in V may now be carried out independently , 
changing the normalization and leaving an integration over only K +1 functional variables 

Zn = C';, yn%a]exp ^^|dr^/Xa(r)p,(r)+ln|||ndr(")exp^-^/Xa(r("))-/3y(r(0),...r(^))j^^ . 

(10) 

The terms Cjv, C^, C'^ are all constants whose value is unimportant. Now in the limit of large N we use that fact 
that 

J drna{r)pa{r) ^ j drMa(r)Pa(r) + 0{-^), (11) 
to eliminate the dependence on the quenched disorder in /9a(r), and obtain 



N 



(12) 

The above integral may be evaluated by the saddle point method, the method works in this limit because the 
thermodynamics is determined solely by the Pa which are the first moments of the random quenched densities Pa- We 
find the following expression for the free energy per number of particles. A'', in each set Sa 



-13 f = max^^ 



j dvY,l^a{v)pa{v) + \n{j dr(")expUEM„(r("))-/3y(r(°),...r(^)) 



+ (13) 



The last term is a beta and independent term and is related to the infinite temperature entropy, it does not 
contribute to the energy of the system, which is what interests us here. The resulting saddle point equation for the 
Ma is 



Pa(ra) = ^ exp(-M„(r„)) / dr(^) exp I - ^ /Xi,(rW) - pV{v 



(0) ... M 



(14) 
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where 



exp 



(15) 



One sees in the above that if /ia(r) is a solution to the saddle point equation then so is ^a{^) + OLa, where Ua 
are arbitrary constants. This is because the constraints Jrfr ^0(1") = 1 are automatically satisfied and so the zero 
frequency Fourier modes of the /i^ are redundant and are zero modes of the theory. We can use one of these free zero 
modes to set Z in Eq. H15() to be equal to one. The saddle point equation then becomes 

Pa(r('^') - exp (-A^alrf"^))) / \{ dr^'^ exp [ - ^ ^^(r^^)) - f3Vir^°\ ■■■ , r^'^), • • • , r^^)) j , (16) 

and we find that the part of the free energy relevant to the calculation of the energy is 

J drJ2Pair)t^a{r), (17) 



and the average energy is given by e = df3f /d(3. This can be manipulated to find that in general, the energy may be 
written: 



e= / Y[dr'^''^Vir'-°\--- ,r'-''\--- ,r(^))exp l~J2fib{r 

•'a \ b 



W^_/3y(r(0),... ,r('^),... ,r(^)) . (18) 



We have verified the formalism for non-uniform distributions in some simple cases, but the rest of the paper will be 
concerned with uniform distributions in domains of unit size. We write Sa(r) = exp(/i£i(r)) to obtain the equations in 
the following form 

«a(r('^)) = JUdr^'^ey^p^-pVir^'K--- ,r('^\--- ,r(^^))l[^-^ 

e = |ndrWv^(r(«),...,rW)n^^oxp(-/?nr(°\---,r(^))). (19) 

a b^a ' ^ ' 

These equations now form the basis for the rest of the paper that explores their consequences. 

In general Eq. pCI) cannot be solved analytically, however in the limit /3 — > cx3 it can be simplified by a saddle point 
method if one writes Sa(r) — exp(— /3K;a(r)) which yields 

?«a(r) = minri...rK.|r„=r{^(ri, • • • ,Ta) - ^Wb{rb)-} (20) 

When we are interested in the maximum rather than the minimal problem we can either change the sign of V or /3, 
changing the sign of (3 simply leads the min above to be changed to a max. In the limit /? — > 00, using Eq. (|17|l . this 
yields the ground state energy or optimal matching cost to be 

^GS^Yl I drwa{r). (21) 

a •' 

The saddle point equations l|19|) are similar to those studied in our work on Hamiltonian paths but with the major 
difference that there are now K + 1 functional order parameters /^a(r). In the case where y is a symmetric function 
and all the Pa are the same for each set, implicitly the case here as we take them all to be uniform, there will clearly 
be a set symmetric solution where /^a(r) — fJ-ij). However it is possible that this symmetry could be spontaneously 
broken. Thus, despite the fact that the matching problem superficially looks somewhat simpler than the Hamiltonian 
path problem, it has the potential to exhibit more complex behavior. 
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The simplest example we can consider is the bipartite matching problem with K = 1. In this case the saddle point 
equations read 



so(r) = J dr 



,ex p (-/3T/(r,rO) 
si(r') 



If one considers the set symmetric solution so(r) = si(r) = stsp(r), the resulting equation is 

J StspyiC ) 

this is exactly the same equation as that occurring for the TSP problem however in the TSP problem stsp(r) = 
exp (/i(r)/2), where /i is again the Lagrange multiplier fixing the density of points and the TSP free energy per site 
is given by 

fiftsp = - j druir) = J dr In (stsp(r)) . (24) 

However we see immediately that for the bipartite matching problem the free energy, within the set symmetric solution 
is given by 



j dv\u{s{v))^(3ftsp. (25) 

Thus the free energy of the TSP with N points is the same as that of a bipartite match between two sets of N 
points. Notice that in the TSP there are N links and in the bipartite match there are also N links, this equivalence 
is conjectured to hold for the two minimal versions in the limit where the energy is scaled so as to become extensive 
as /3 ^ oo. If the set symmetric solution is indeed that valid everywhere, our result shows the strict thermodynamic 
equivalence of the two models at any finite temperature and for an arbitrary distance function ! 

For pairwise potentials arranged as paths Eq. Q and uniform probability density over the domain we may search 
for general K solutions along the lines of the one found above for the bipartite case. The non-cyclic case turns out 
to be more tractable. We search for a solution in which Sa{v) — exp(/ia(r)) is the same for all sets except the first 
(a = 0) and last (a = K) where the ends of the matching path lie. 

r^(r) « = l-^-l (26) 
^ ' \ se{v) a = 0,K ^ ' 

The saddle point equations now read 

r exp -/3Eti^p( 

SE{ro) ^ dri...drK ^- . . / (27) 

J s{ri)...s(rK-i)sE{rK) 

r exp (-/3X;f=i Vp(ra-i,ra)) 

s(ri) = / drodrs . . . dvK \ , , , T^^' (28) 

J SE{ro)s(r2) . . .SEirK) 

and free energy corresponding to this saddle point is 



f3f 



j dr [2 \n{sE) + {K-l) ln(s)] . (29) 



Using the structure of the pairwise potential, these equations can be solved in terms of a set of K coupled integral 
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equations for quantities ia(r) 

ti{r 

tK-2{r 



dr' exp {—f3Vp{r, r' 
dr'exp(-/3T/p(r,r' 
dr'exp {-pVp{r,r' 
dr' exp (— /3Vp(r, r' 
dr'exp {^pVp{r,r' 
dr'exp(-/3V^p(r,r' 



tK-iir') ^tK-i 

tK-2{r') Hk-2 



ii(r')iK-2(r') 
tojr') 

t2{v')tK~2{v') 

^o(r') 



^0 



to 



t2tK-2 ' 



tK -sir' )tK-2ir') 
1 _ . 1 ^ 



to 



tK-3tK-2 



) 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



where we have introduced T as the integral operator appearing throughout. Once the taS are determined the solution 
is given by: 



S£;(r) = tK-i{r) 
s(r) = to{r)tK-2{r). 

This set of equations always has a solution with all ia's taken the same ta(r) = stsp(r), this then yields 

ssir) = stspir) 
s(r) = sLp(r), 

and inserting this into the free energy we have 

(3f = -2K [ drlnistsp) - KPftsp. 



(36) 
(37) 



(38) 
(39) 



(40) 



The interpretation is clearly that the free energy per link of this matching problem is exactly the same as for the TSP. 
This result is supported by direct Monte Carlo simulation in the case of a two dimensional box. This is a generalization 
of the thermodynamic equivalence seen between the TSP and the bipartite matching problem. Unfortunately, neither 
for closed paths nor for fully connected matches have we found any general results and we shall consider only special 
cases in the following section. Finally we comment that for sufficiently symmetric potentials on closed domains (e.g. 
with periodic boundary solutions) we find the solution Sa = Constant for all a. This means that the annealed 
approximation to free energy is exact and it is easy to show that ecs is just the ground state energy of a single 
molecule or match with m = K + 1 sites in the case where each of the m sites is free to move. 



III. TRIPARTITE MATCHING- NUMERICAL RESULTS 



Our analytic efforts have not revealed a solution for a symmetric K = 2 tripartite system other than a constant on 
symmetric closed domains. We have therefore performed some numerical simulations of the triangle potential. We 
have considered two domains each with boundaries: a one dimension line and in a 2D box. After preliminary tests to 
estimate the strength of 1 /N corrections and the time necessary to equilibrate we choose the following Monte Carlo 
parameters: N = 1000, equilibration for one million steps and measurements taken over the succeeding one million 
steps. Energies are averaged over 100 samples with different sets of random points. Our codes are checked by showing 
that the bipartite K — 1 case agrees with the TSP result. 

In Fig. Q the Monte Carlo results for the energy are confronted with a numerical (iterative) solution of the 
equations. We find that in all cases, the iterative approach converges to a symmetric solution with all Sa (r)'s the 
same, even when each Sa(r) is seeded with rather different initial conditions. In two dimensions, the iterative solution 
is slow for a reasonable discretization of the domain, however we find that for the temperatures we consider the 
solutions are smooth and the procedure gives good results even for a discretization with only 400 points. 

We have also considered tripartite matching with the star potential and again see no evidence for breaking of the 
symmetry between the sets. 
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FIG. 5: Tripartite matching with a triangular potential in 2 dimensions. Points are from Monte Carlo (the error bars are too 
small to show) and the continuous line is an iterative solution of Ea. Negative temperatures correspond to the maximal 
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FIG. 6: The function 'w{x) for tripartite matching on a one dimensional line. Obtained by numerically solving Eq. 1)19^ for 
/3 = 70.0. The parameters agree with the solution based on a greedy heuristic discussed in the text. 



IV. ZERO TEMPERATURE RESULTS IN ONE AND TWO DIMENSIONS AND CORRESPONDING 

HEURISTICS 

In this section we consider one and two dimensional examples where our formalism can solve the original maximal 
optimization problem. We study Eqs. H2U|I . (|21l) and find that, as was the case for the TSP i9|], the method of solution 
throws light on local heuristics which might be used to solve specific instances. 

First consider maximal tripartite matching in 1 dimension with a triangular potential based on a pairwise potential 
that is simply the distance between the points. 

Vtriangieixo,xi,X2) = l^^o - Xi\ + \xi ~ + \x2 - | = 2 [max(xo , Xi , X2) - min(xo , Xi , X2 )] . (41) 

Evidently the location of the middle point does not contribute to the potential. We search for a piecewise linear 
function w(x) consisting of three linear sections evenly dividing the unit interval as shown in Fig. I^. A l/3rd jump 
heuristic guesses that the maximum occurs for points arranged as shown in Fig. iQ. This basically means that the 



X 


x+1/3 x+2/3 




1/3* i/3* 1 



FIG. 7: Arrangement of points for the one third jump heuristic. 
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X 2/3-x 1-x 

kr* — 1/3 • i/3 • 

FIG. 8: Arrangement of points for the greedy heuristic. 



match consists of a central point in the central third of the interval, matched with points in the first and last third, 
each at a distance of a third from the central point. In this case the Eq. (|20|l can be written: 

W{X)+W{X + 1/S)+W{X + 2/3) ^ Vt„angle (x, x+l/3,x + 2/3) = 4/3, (42) 

where a; : < a; < 1/3. This of course immediately gives the average value of 4/3 per match in the optimal match. 
We now make a piecewise linear anzatz for w{x) (which is symmetric about a; = i) 

w{x) = a\x-^\+b forxe [0,i]U[^,l] 

= a'\x-^\+b' forxe[i,^]. (43) 

Immediately one sees from Eq. (|42ll that we mus have a' — 0. Now substituting this into the Eq. (|20|l we find 

^a + 2b + b' = ^, (44) 

The coefficient a can be determined by looking at locally at Eq. ipn|l when a; G [5, |] and writing xi = a; — ^ + ei and 

X2 = a; + -i + £2 5 this gives 



b' = max I ^ + 2(e2 - ei) - a{e2 - ci) - - 2b\ 

«1,«2 {6 3 J 



(45) 



and the local stationarity with respect to ei and £2 then gives us that a = 2. Now we assume continuity of w which 
gives 

b' = l (46) 

This thus yields the solution for Eq. (|20|l in complete agreement with the numerical resolution of the low temperature 
saddle point equations shown in Fig. ©. 

An alternative heuristic for this triangular distance potential is the greedy heuristic with the arrangement of points 
shown in Fig. ijSJ. In this case the maximization equation is: 

w{x) + w{2/5 ~x) + w{l -x) = Vt„angieix, 2/3 - X, 1 - a;) = 2(1 - 2a;). (47) 

Interestingly despite the fact that the heuristic here is different it yields the same function w{x), and consequently 
the same average value per match of the optimal match. This is presumably a peculiarity of linear pairwise potentials 
in one dimension and a similar phenomenon is seen for the Maximal TSP j9(. 

These issues generalize to other cyclic (and possibly fully connected) potentials in 1-dimensional m-partite matching. 
The star potential based on the same distance metric, does not lead to such a simple solution as the locations of all 
three points are important and w{x) is no longer piecewise linear. 

Now, remaining in 1-dimension, consider the triangle potential built from a quadratic pairwise potential. The 
location of all three points are needed to determine the cost of a match (in contrast to the linear distance function) , 
and it is harder to guess a heuristic. However, numerical experiments hint that the solution w{x) is quadratic so we 
try the anzatz w(a;) = a + b{x— i)^ (which respects the symmetry about the mid-point x = 1/2) in the maximization 
equation: 

w{x) — max \{xi — a;)^ + (3:2 — a;)^ + (xi — ^2)^ — w{xi) ~ ^(3:2)] • (48) 
The maximization equations for right hand side yield 

(2 - b)xi - X2 = 

{2-b)x2-xi = x-^b. (49) 



A manifestly non-optimal solution to the above is X2 = and to avoid this we require the above equation to have 
more than one solution, this means that 5 = 1 or 6 = 3. The choice 6=1 yields x = 1/2 which is clearly generally 
not the case. The choice 6 = 3 yields x + xi + X2 = S/2 which can be written as 

^{X + Xi + X2) = ^. (50) 

meaning that the center of mass of optimal triangles is in the center of the interval. Using this condition in the 
maximal equation, all parameters are determined, we find that a = and the energy is given by 3/4 which can 
further be checked by numerical iteration of Eq. 116|l . The condition above is not a full heuristic in that it does 
not determine the location of all points in the match when given just one. However, it is supported by Monte Carlo 
simulations where the sum of the coordinates of all matching triangles is a bell shaped distribution with correct mean, 
and variance decreasing rapidly with temperature. We expect this result to generalize to m-partite matching with a 
cyclic potential with quadratic pairwise Vp. It also generalizes to a star potential with quadratic Vp where we find 
energy 1/4. Moreover, although it is not an exact solution, the quadratic form of w{x) is a good approximation even 
a quite elevated temperatures. 

The amenability of quadratic potentials carries over to two dimensions. We consider tripartite matching with a 
triangle potential in a disc so as to preserve rotational symmetry and take w{r) to be independent of angle. The 
function to be maximized is smooth and by differentiating we find that a symmetric maximum occurs when the center 
of mass condition ri -|-r2-|-r3 = holds. As in the 1-dimensional case this is not a full heuristic, however it is sufficient 
to determine the energy which is 9/2tt. This is consistent with numerical iteration of Eq. (|16|l at low temperature, 
though in 2 dimensions numerical resolution of the saddle point at low temperatures is difficult. Furthermore Monte 
Carlo simulations indicate that the center of mass condition indeed becomes sharper as the temperature is reduced. 
A specific heuristic which yields this energy makes greedy matches based on equilateral triangles with center at the 
origin. The same heuristic can be used in the case of a triangle potential in 2-dimensions based on a linear distance 
pairwise potential; it leads to an energy of 2^/3/t:, but although this is in good agreement numerical results we have 
not been able to prove hat it is the true maximum. 

V. CONCLUSIONS 

We have presented a statistical mechanics based approach to multi-index matching problems with site disorder, 
specifically when the points to be matched are randomly distributed in a domain V. The cost functions considered 
are functions of the relative distances between the points and are consequently correlated, in contrast to the random 
link model. We have analyzed the m-partite version of these problems where each match contains one element from 
m distinct sets, each of which can in principle have different probability distributions in the domain V. In the 
thermodynamic limit the system is described by m functional order parameters. We have concentrated on maximal 
matching problems as the ground state energy in this case is extensive with the scaling we employ. 

In the special cases where the cost function is sufficiently symmetric and the probability distributions of each set 
are the same we find set symmetric solutions where the functional order parameter of each set is the same. We have 
not found any evidence for a breaking of this symmetry, an open question is to whether there are geometries (e.g. 
effects of boundaries and spatial dimensionality) where this symmetry is spontaneously broken. 

For open chain potentials (which of course includes bipartite matching) we find a solution to the saddle point 
equations which gives exactly the same free energy as the traveling salesman problem with the same number of 
links. Furthermore this equivalence holds at all temperatures and is independent of the precise functional form of 
the pairwise potential from which the chain potential is constructed. Interestingly this observation is analogous to an 
equivalence which is conjectured to hold between Euclidean bipartite matching and the Euclidean TSP in the zero 
temperature limit of the minimal problem where the cost function is scaled to give an extensive ground state energy. 

As in the TSP we find that the zero temperature saddle point equations can be solved via ansatze inspired by 
heuristics to find the optimal match. We have been able to analytically solve the saddle point equations for the 
function u; in a number of instances of maximal matching in one and two dimensions. Also we have numerically 
verified that the average optimal match is accurately predicted and that the heuristic which was used to solve the 
saddle point equation is indeed that associated with the optimal match. Finally we remark that knowledge of the 
function w does not completely specify a heuristic however it does give some partial, potentially useful, information 
about the heuristic. Indeed we have demonstrated that two completely different heuristics can give the same average 
of the optimal match for a simple Euclidean tripartite-matching problem in one dimension. 
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